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Abstract 

iy~j ' A study of the high angular momentum particles 'atmosphere' near the Schwarzschild 

, black hole horizon suggested that strong gravitational interactions occur at invariant distance 

of the order of \/M. We present a generalization of this result to the Kerr-Newman black 
hole case. It is shown that the larger charge and angular momentum black hole bears, the 

qq ' larger invariant distance at which strong gravitational interactions occur becomes. This 

■ invariant distance is of order ^/r^/(r+ — r_). This implies, that the Planckian structure of 

O the Hawking radiation of extreme black holes is completely broken. 

0" 1 

U 



PACS: 04.60.-m; 04.70Dy 

Keywords: black hole, atmosphere, fluctuations, semi-classical approximation 



X 

1 Introduction 



The black hole radiation is a direct consequence of non-static character of the collapsing body 
metric. A wave propagated through the collapsing body becomes red-shifted, hence an in- 
vacuum state defined with respect to modes at past null infinity differs from an out-vacuum 
state, defined with respect to modes at future null infinity. The red-shift is exponentially large 
if collapse leads to the black hole formation. 

It was shown by Hawking^ that the black hole radiation spectrum is the Planckian one, 
provided that there is no back reaction of emerging quanta on metric; the matter fields propagate 
on a classical background metric, and the wave equation is valid at all scales. 

This semi-classical approach of Hawking becomes doubtful at scales of the order of l p = 
\J Gh/c 6 , i.e. the ultraviolet part u> > lp 1 of the Hawking radiation spectrum is possibly non- 
Planckian. However, such huge energies are of little interest from an experimentalist point of 
view. It turns out that problems begin at much smaller scale. 
1 e-mail:tuchin@post. tau.ac.il 



1 



We can gain more by putting in question the possibility to neglect the radiation back reaction. 
Existence of energy flux from the black hole at infinity implies that the black hole mass decreases. 
However, as long as the mass of the black hole is large compared to the Planck mass, the 
rate of evolution of the black hole is small compared to the characteristic time for the light 
to cross the gravitational radius ||. Thus, we can describe the black hole by a sequence of 
stationary solutions; in each solution the back reaction influence on the radiation spectrum can 
be neglected ]|]. 

On the other hand, despite the fact that the Riemann tensor is small near the black hole 
horizon R ~ 1/M 2 the normal modes get exponentially red-shifted^]. This implies, that strong 
gravitational interactions occur there and, possibly, alter the Planckian character of the spec- 
trum. 

The brick wall model proposed by 't Hooft is a simple model that provides some insight 
into the problemQ. The idea is to calculate the number of the (scalar) wave equation solutions 
near the horizon of the Schwarzschild black hole. This number can be interpreted as a partition 
function of the black hole " atmosphere" — a gas of high angular momentum Hawking particles 
reflected back by the space-time curvature[Q]. The atmosphere's entropy turns out to be diver- 
gent at the horizon. However, one can introduce an invariant cut-off p of the order of 1 and 
independent of M, such that the atmosphere's entropy becomes exactly the Bekenstein-Hawking 
one. This implies that the strong gravitational interactions occur at invariant distance of the 
order of 1. Several other authors arrive at this conclusion using other models ||. 

However, some authors arrive at opinion that p is greater that 1@. In particular, Casher 
et. al. showed that p is of the order Af 1//3 , provided that one takes account of the atmosphere's 
thermal fluctuations. This implies that the Hawking radiation Planckian structure is broken 
at uj > Af -1 / 3 . Therefore, unlike the Planckian spectrum quanta, the real black hole radiation 
quanta are correlated, and information about the state of in-falling matter gets encoded out of 
the horizon J2|. 

The brick wall model may be criticized for it deals with non-renormalized stress tensor. 
We can avoid this difficulty as discussed in ref. ||. Both approaches yield the same result for 
the atmosphere's entropy. So, we shall apply the brick wall model again in order to obtain an 
expression for the charged and rotating black hole entropy (sec. [| and 2.3)|11|. 

Along the way we shall learn more about the atmosphere. In particular, it will be shown 
that the main contribution to the partition function comes from the particles emitted with a 
charge and angular momentum of the same sign as the black hole ones. 

In section 2fl the atmosphere thermal fluctuations will be studied. It will be shown that 
strong gravitational interactions near the charge black hole horizon occur at an invariant distance 
of the order of [r+/(r + — r-)] 1 / 3 . In section 2J3 it will be argued that this result is valid for the 
Kerr black hole also. 



In section g we shall confirm result of section by applying the shock wave model [12] to 
study the gravitational interactions between in- falling and out-going particles. 
We discuss the results in section ^ 
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2 Thermal Properties of Atmosphere 
2.1 Charged Black Holes 

The charged black hole metric (Reissner-Nordstr0m geometry) is given by: 

r z D 

where D = r 2 — 2Mr + Q 2 , M and Q are the mass and the charge of black hole respectively. 
The event horizon is located at surface r + defined as 



r+ = M+^M 2 -Q 2 . (2) 

A distant inertial observer observes the Hawking radiation — energy flux from the black hole. 
This flux has a Planckian structure provided that we neglect the back reaction of this radiation 
on the metric]!]]. Therefore, the Hawking radiation can be considered as a black body radiation 
in a thermal bath with the temperature Th 



where A = Anr 2 ^ is the horizon area|8|. 

Consider a massless neutral scalar field 4> i n the charged black hole background. One can 
separate variables in the wave equation 

g^v^A = o 

and obtain the following expressions for the basis functions ||: 

UEi m (x)=r- 1 R l (r)Y lm {e,i P )e- iEt , 

where I = 0, 1, 2 . . ., m = —I, — I + 1, . . . , I — 1, 1, Y[ m (9, ip) are the spherical harmonics, E is the 
hamiltonian eigenvalue. The radial part Ri{r) satisfies the following equation: 

[^-V El (r))R l (r) = 0, (4) 
where r* and Vei(t) are defined by the following relations: 



and 



D 2MrD 



V El {r) = -E 2 + l 2 - i + —^. (6) 
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The centrifugal barrier is 

It is attractive for r + < r < r% and repulsive for r > n, where r\ is the root of the equation 
Wl -o- 



n = !Uf + ^-!<yU, (7) 

The tunneling through the angular momentum barrier may be neglected for all but the lowest 
angular momentum modes since Vf ~ ^ 2 0- 

We see that near the horizon there exists an atmosphere of the neutral high angular momen- 
tum particles. In general, the same is also true for the charged scalar particles which are defined 
with respect to basis functions of the following equation: 

g^(V^-ieA^(V v -ieA v )(l>. (8) 

Let us study properties of this thermodynamical system. At first we shall calculate the number 
of the wave equation solutions T. As long as M 3> 1 we can rely on WKB approximation. So 
we get: 

r = (2vr)- 3 £ f dx k d Pk ; k = r,9, <p, (9) 

9 J 

where the sum runs over all degenerate energy levels. 

The covariant radial momentum p r associated with the differential operator p r in the wave 
equation is given by|| [7j]: 

p 2 r = D- 2 (r 2 E-eQr) 2 -D'H 2 . (10) 

Here e is the charge of the scalar particles. Since the metric is spherically symmetric integration 
over the 6 degree of freedom in eq. (g) is trivial. Thus, 



T(E) = ^2(2tt)- 2 J2f dl ff dr f d P r ~^J dl ( 21 + 1) J dr J dE 



dp r 

1e' 



This expression is divergent at the horizon and at infinity. To proceed further we should intro- 
duce cut-offs. To this end let us use the brick wall model H. 

We introduce the following boundary conditions: 4>(x) = if r > L, and <fi(x) = if 
r < r+ + h, where L and h are the infrared and ultraviolet cutoffs respectively. We are interested 
in the contribution of the horizon, i.e. in the term 0(/i _1 ). The term proportional to L 3 is the 
usual contribution from the vacuum surrounding the system at large distances and is of little 
relevance here. So we find: 

D ss (r+ - r_)/i, (11) 
1 r ftm ri(E - e$)D- 2 h 



l r f m 
T{E) « - dE dll 
vr J Jo 



D- 2 r%{E - e<&) 2 - D-H 2 
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dE(E - e$) 2 , 



irh (r+ — r_) 2 

where l m is the largest possible value of I which arises from the obvious condition p 2 > 0: 

r\(E-e<$>)^ 
~ h(r + — r_) 2 ' 

and $ = Q/r + is the electric potential at the horizon. 
The free energy of the atmosphere is given by: 



(12) 



F = T H J dT{E) ln(l-e (e *- £)/3 ). 



(13) 



This integral exists only if the condition e<l? — E < holds, i.e. the supperradiant modes do not 
contribute to the atmosphere's free energy. Hence, 

1 r 6 
T(E) = — - x 



irh (r+ — r_) 2 

rE 



X 



i-E eE 

\ (E' - e&) 2 dE'6{-eQ) + / {E' - e<&) 2 dE'6{eQ) 

JO Je$ 



and also 



F = 
+ 



3irh (r+ — r_) s 



{(£ - e$) 3 + (e$) 3 }6(-eQ) + {(E 1 - e$>f}9{eQ) 



= T(E)_+T(E) 



+ 



T H T(E) ln(l - e^- E ^)\% =0 - / dE 

Jo 

T H T(E) ln(l - e^- E ^)\% =e „ - / 

-l^ dE ^Shi e{ - eQ) -L 

F- + F+, 



T(E) 



6(-eQ) 
9{eQ) 



e $ e (s-e*)/3 _ x 



where we have assigned the subscript + to the contribution of particles with the same sign of 
charge as the black hole's one, and the subscript — to the contribution of their oppositely charged 
antiparticles. The step-function 6(x) is defined as usual: 9{x) = 1, if x > 0, and 6(x) = 0, if 
x < 0. Substitution of F yields: 

F_ = 



3irh(r + — r_) 2 Jo 



e (E-e<S>)f3 _ 1 



F, = -- 



+ 



L 



(E — e&) 3 
dE m „ ) r , eQ > 0. 



3nh(r + - r_) 2 J e $ e^-^ - 1 



(14) 
(15) 
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As long as we assume that M>1, condition /3 3> 1 holds for the black hole of any charge. Then 
neglecting — 1 in the integrand denominator and changing the variable of integration E — e<£ = y 
one gets: 



r 6 



= - 3jft(r f_ r ,2 e,StT H - teP* + 6) . e<3<0; (16) 

F + = ~TTT^ ^ / V^~ vPd V = "TTT^ ^ 6T "' e ^ ^ a ( 17 ) 

37r/i(r + — r^y Jo 6TTh[r + — r_y 

From eqs.(16| 1?]) we deduce that the contribution of particles carrying a charge with sign opposite 
to the black hole's one is negligible compared to the contribution of their antiparticles. This 
is because they are pulled into the black hole by the electrostatic field as soon as they emerge 
and thus spend only a short time outside the horizon. Conversely, particles with the same sign 
of charge as the black hole's one are pushed out of the horizon, then scatter off the effective 
potential V\ (see eq.(^)) and then the lowest angular momentum ones escape to infinity, and the 
others return to the atmosphere. 

Note that F + does not depend on e. Thus, the thermal properties of the charged atmosphere 
(and, hence, contribution of the horizon to these properties) are the same as the neutral one 
and defined by the black hole mass and charge completely. 

In the limit of the neutral scalar field (e = 0) eqs. (17,16) reduce to 



2r° 

irh(r + — r_) 2 



F n = F_(e = 0) = F + = - — + _ 1%. (18) 



The entropy of the neutral atmosphere is given by: 



_ dF n 8(r + -r_) 



This coincides with the Hawking-Bekenstein entropy S = \A (T|, |5| for the following value of the 
cutoff: 

8(r + — r_) 

(4 7r )3 7r 2 r 2 



h = K= (20) 



The fact that h depends upon M and Q is merely due to the special choice of coordinates. 
Define the invariant distance as follows: 



P 



r=r + +h rr=r + +h 2r i i— 

ds= yj-g rr dr fx - TTJ^vh- ( 21 ) 



So, Pn = P {h n ) = ir-V 2 !- 1 / 2 . 

We see that p is a property of the horizon independent of M and Q. The same result for the 
Schwarzschild black hole was obtained by 't Hooftp]]. 
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In terms of the invariant distance the neutral atmosphere entropy reads: 



32 r% 
(4vr) 3 vr p 2 n 



\eQ\P ~ \eQ\- — ~ \eQ\ M — == » 1; 



Suppose now that e / 0. In any realistic black hole the following inequality holds: \eQ\/3 S> 1. 
Indeed, 

M 2 . n . M 

~ \eO\ — - 

'(r+-r_) 1 ^ Vl " (<9/M) 

which is clearly true since Q 3> e e i ec t ron ~ 10 _1 and M>1. Hence, in the leading order eq. (p^) 
gives: 

F- = ~ , ( r&+ . a 2f^(eS) a . 
irn{r + — r_J z 

Since -F + does not depend on e the entropy of the charged atmosphere has only an exponen- 
tially small dependence on e: 



S c — *S*-)- -\- S— 

32 ri ri |e$| 



(4-7r) 4 7r p 2 47r/i (r + — r_) : 



e- |e * l/3 « 5 n . (22) 



It is seen that at p c = p n this expression returns to the Hawking-Bekenstein formulae, as it has 
to, because, as we see, the thermal properties of the atmosphere, which were used to define the 
cutoff, do not depend on the scalar field charge e. It will be shown in the next section that the 
invariant distance in the Kerr metric is also independent of the atmosphere particle's angular 
momentum projection on the symmetry axis m. 

For non-extreme black holes of mass greater than 10 15 g it is a Klein-paradox process that 
dominates the charged pair production, and the emission rate is governed by a Schwinger-type 
formula (for pair production in a constant electric field). For black holes of smaller mass Hawking 
thermal process dominates §. 

2.2 Horizon Fluctuations 

It was suggested in M that due to the atmosphere of high angular momentum particles, strong 
gravitational interactions occur near the Schwarzschild black hole horizon at an invariant dis- 
tance of the order of Af 1 / 3 . This is because the total energy of the black hole's atmosphere 
fluctuates as any thermodynamical variable. Let us extend this idea to the case of the charged 
black holes. 

According to the rules of thermodynamics, given the energy of the atmosphere U (p) between 
the surfaces r = r + + h{p) and r = r%, and the number of particles N(p) between these surfaces, 
the black hole's mass fluctuation between r = and r = r + + h(p) is 

AM(p) = AU(p) ~ -^L=, (23) 
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provided that Riemann tensor is sufficiently small near the horizon R ~ M~ 2 <C 1, so that we 
can apply the usual flat spacetime thermodynamic rules in the vicinity of the horizon. We have 
used the fact that the total energy M defined in (jl|) is fixed from the point of view of an external 
observer. The energy of the atmosphere above the surface r = r + + h(p) is given by: 

The number of particles is N ~ S. Therefore: 

AM( P )~ (r+ - r - ] 



r+P 

Right now we are facing problem mentioned in the section EL we had to use renormalized 
total energy JJ ren instead of U . However, as far as only the variation of the total energy is 
needed, one may replace AU ren by AU. Indeed, all divergent terms in the energy-momentum 
tensor are functions of the curvature tensor components and their derivatives only|Q, so they 
are canceled from the expression for fluctuation of the total matter energy^. 

Let us estimate how the fluctuating mass gives rise to uncertainty in the location of the 
horizon. A point r' is outside the horizon if r' — r + (M(r')) > 0, where M(r') is the black hole 
energy between r = and r = r' . Clearly, if 

A(r' - r+(M(r'))) = A(r+(M(r')) > r' - r+(M(r')) = h, (24) 

then the point r' is in the superposition of being inside and outside the horizon. Using eq. (||) 
we obtain: 

2AM r + 1 , x 

Ar+ = = -. (25 

r + — r_ p 



From eq. (21) we see that relation Ar+ < h holds if 



p r\ 



that is, the semi-classical approach is valid as long as 



.,2 



1/3 



P > Pmin = • (26) 

V r + — r_ J 

This equation implies that in Schwarzschild metric p m in = M -1 / 3 , see ref. Q. In the case of 
extreme black hole p m in becomes infinite. 



2 In the ref. J^j] S.Mukohyama and W. Israel argued that if one correctly identifies the ground state — the 
Boulware vacuum, then the total energy near the black hole horizon is finite. They use this fact to advocate 
the brick-wall model. Conclusions of this section are independent of choice of the ground state; the thermal 
fluctuations depend only on the number of modes near the horizon. 
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We have derived eq. (25) neglecting fluctuations of the atmosphere's total charge Q. Never- 
theless, it remains valid if take account of the Q-fluctuation. Indeed, 



Ar + - \ / ( — — — ) 2 AM 2 + - — — — ^2 AQ 2 



= ac/ 2 + ?aq 2 = u-^q . 

r + — r_ y rfj_ r + — r_ V r + J 

Since the quantity in the curly brackets of the last equation is additive we estimate it's fluctuation 
as follows: 



V r + J r + p 



Putting all this together leads to eq. (| 
2.3 Kerr Black Holes 

Consider a rotating black hole with projection of angular momentum on the symmetry axis J 
(called simply black hole's "angular momentum"). The metric in Boyer-Lindquist coordinates 
is given byQ: 

ds 2 = (l- dt 2 - ^dr 2 - £ W 



£ 2 J D 

2Mra 2 . „\ . Q „ , , 2 • 2Mra 
l] 2 



r 2 + a 2 + ^ 2 sin 2 9 sin 2 9d<p 2 + sin 2 0d<pdt. (27) 



where a = J/M, 

D = r 2 — 2Mr + a 2 , S 2 = r 2 + a 2 cos 2 0. (28) 
The event horizon r + and surface r_ are defined as 

r+ = M+VM 2 -a 2 . (29) 

There is a surface such that no static observer can exists inside it (frame dragging effect). 
Definition of this surface, called static limit, is 



r (9) = M + VM 2 - a 2 cos 2 9. 

As in the case of the charged black hole there exists Hawking radiation with the temperature 
T H given by 0: 

_ (r + - r_) _ (r+ - r_) _ (r + - r_) 
H ^ 4vr(r^ + a 2 ) 8vrAfr + 

Inertial observer at infinity measures flux of energy and angular momentum from the rotating 
black hole. We are interested to study what happens to the system of Hawking particles reflected 
back to the black hole by the space-time curvature. 
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Recall that in the Kerr metric one can separate variables in the wave equation (||)|| |i~0[ . 
The basis functions are: 

u E x lm m(x) = (r 2 + a 2 r 1 / 2 R lm (r)S lm (co S 0)e im ^e^ Et , 

where Si m is the spherical harmonics with eigenvalue Xi m (aE), I = 0,1,2..., m = + 
1,... 1-1,1. 

The radial part Ri m (r) satisfies the following equation: 

d 2 

V E im{r)}Rim{r) = 0, 



(30) 



dr* 2 

where r* and Ve/ m (r) are defined by the following relations: 

dr* r 2 + a 2 
~dr~ = D ' 

V E im(r) = -(E-m a ) 2 + \i m (aE) D + 

r A _|_ a z _|_ a zy 

2{Mr-a 2 )D 3a 2 D 2 
(r 2 + a 2 ) 3 + (r 2 + a 2 ) 4 ' ( } 



The centrifugal barrier is 



Vkmir) = km{aE)- ° 



(r 2 + a 2 ) 2 

It is attractive for r + < r < r\ and repulsive for r > n, where ri is a physically acceptable root 

dr 

n = M + 2JM 2 --a 2 -cos-, (32) 
V 3 3 



of equation = 0: 



where 

1 - a 2 /M 2 



COS Q 



(l-a 2 /3M 2 ) 3 /2- 

Once again the tunneling through the angular momentum barrier may be neglected for all but 
the lowest angular momentum modes since Vf m ~ A; m ~ I 2 . 

The thermal atmosphere of the black hole extends from the horizon up to the surface r = r\ , 
where r\ is the turning point of the centrifugal barrier. Study of two limit cases a = and 
a = M shows that r\ > ro for any black hole. In other words the ergosphere lies under the 
external boundary of the atmosphere. 

Let us calculate a number of wave equation solutions in the vicinity of the black hole. This 
number is equal to the twice of number of quantum states of the atmosphere, because of two-fold 
degeneration of energy levels (along the symmetry axis). In the WKB approximation: 



T = 2(27t)~ 3 ^ J dip J dr J dp r J d9 J dp e . 
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Here the covariant radial momentum p r associated with the differential operator p r in the wave 
equation is given by H: 



P 



2 



D- 2 [(r 2 + a 2 )E-am} 2 - Xi m D-\ (33) 



and the polar momentum pg associated with pg: 



2 



\ lm -(aE sm0-— -f. (34) 
smt) 



Introduce a cut-offs h and L which are defined as <j)(x) = if r > L and r < r + + h(0). Now 
allow /i to be some function of 6. 

Let us change variables of integration replacing the pair (p r ,pe) by the pair of constants of 
motion (X,E) (here A; m = A ). Jacobian of this transformation is 

d(p r ,Pe) 



d(X,E) 

2D~ 2 {r 2 + + a 2 ) 2 (E - mQ H ) - 2D- 1 (aEsm6 - m/sm6)asm6 
4 Pr (\,E)pg(\,E) 

Here £Ih = a/2Mr + is the angular velocity of the black holeQ. Since the main contribution 
to the atmosphere's partition function arises from the region near the horizon (r ~ r + + h) we 
neglect the second term in the last equation and obtain (note eq. (|i"l|)) 

,d(pr,Pe) 



m 2(27r)- 3 2-KhJ2 J d6 J dX J dE x 



D- 2 {r\ + a 2 ) 2 (E - mQ H ) + O^- 1 ) 



2J2D- 2 {r\ + a 2 ) 2 (E - mO ff ) 2 - XD^y/X - (aEsm0 - m/sin6<) 2 



where D near the horizon is given by ( JX1| 
Now, the free energy for large (3 is 



F = T H Y,J dT(E,m)ln(l 



e {mU H -E)f3 



h T H {r 2 + + a 2 ) 2 f f^ m [°° xe-' xl3 
da I aX I ax- 



(2vr) 2 D J J ' Jo ^( r 2 + a 2)2 x 2 _ XD 

x [dm — 1 + o( e -M"«/?) 

7 — [asm6(x + ttiCIh) — m/sin8] 2 

-^T H {rl + a 2 f f d6sm9 H dxx 2 e~ xl3 D~ 2 

2lT J Jo 

1 ^ 4 id + a 2 ) 3 /■ sin0 
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We have introduced a new variable of integration x = E — mQ,H and taken account of the fact 
that contribution of particles with m < is of the order of 0(e~\ m ^ Hl3 ) (like in the charged 
black hole case). The radial momentum p r vanishes at A = A m . 
It is easy to convince, that the entropy is given by 

s.iu&i+f^r**?. (35 ) 

7r r + -r_HJn h 



This formula is consistent with eq. (|l9|). (To see this put Q = in eq. (19) and a = in 
eq. (^)and compare). 

Let us define invariant distance p as in eq. (|2l|): 



ir\ + a 2 cos 2 i 
r + — r_ 



h{6). (36) 



Substitution to (p5|) gives 



5 = 71 Ns ( r +~ r ~) 



d#sin6< if 1 d£(rl+a 2 ^ 



7r(47r) 3 7o ft 7r(47r) 3 7-1 /) 2 (£) 

We may choose h(8) in such a way that this equation will coincide with the Hawking- Bekenstein 
expression for the black-hole entropy. For example, take h(M,a\6) = h(M, a)/(cos~ 2 8 + 3). 
The problem is that there are infinitely many ways to do this, and thus, one cannot fix the 
6 dependence of h. Despite this difficulty, however, we see that whatever function h(9) that 
matches the Hawking-Bekenstein expression one chooses, the invariant cutoff p does not depend 
on the black hole parameters M and J and on the atmosphere's thermal properties. It is 
characterized by the event horizon only. 

The fact that we do not know the exact form of the function h(8) does not matter, if we are 
interested only in ratios of such thermodynamical variables as the total energy, entropy etc. In 
this case we need not know this function. This is the case if we are going to calculate the horizon 



fluctuations. We repeat the by now familiar procedure of sec. 2J2 and arrive at the formula (|2q) 
with r+ and r_ defined by eq. (|29|). Note that deriving this result we have effectively averaged 
thermodynamical variables over all possible choices of h which yield the Bekenstein-Hawking 
entropy. 

It is straightforward to generalize this result to the Kerr-Newman black hole. We will refer 
to eq. (p6|) as such generalization with 



r+ = M+^M 2 -a 2 -Q 2 . (37) 

It should be emphasized, that despite the fact that the p m in becomes infinite for the extreme 
black holes, the condition h <S M under which all calculations were done still holds. Indeed, 
using eq. ( p6| ) this condition reads 

P « (M 2 -a 2 -Q 2 )V4- ( 38 ) 



It is seen by substituting p m in from eq. (p6[) into eq. (|3g|) that the p m in always matches (38) 
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3 Shock Wave Model 



In this section we shall show that strong gravitational interaction in the vicinity of the charged 
black hole occurs at p given by eq. fl2q) by studying the gravitational interactions between an 
in- falling particle and the thermal atmosphere. 

It is shown in the Appendix, that propagation of a massless particle on the charged black 
hole background can be described by the shock wave of a special form, so that the resulting 
metric is given by eqs. (|48|),(|57D,(|58|) and (|60|). Consider some massless particle falling into the 
black hole. The atmosphere (Hawking) particles interact with the shock wave generated by 
this particle. What is the probability that the state of the atmosphere is the same after the 
interaction took place? 

At first, we shall study interaction of the in- falling particle with the one test atmosphere 
particle. Note that the Reissner-Nordstr0m geometry (Q) takes the same approximate form as 
the Schwarzschild one near the horizon: 

ds 2 = p2{r+ ~ r - )2 dt 2 - dp 2 - dx 2 - dy 2 , 
4ri 

provided that we study a region with transverse distances much smaller than M . Here p is 
defined in (^) independently of the black hole charge. Denote 



M 



2(r+ — r_) 
then this metric reads 

ds 2 = -^dt 2 - dp 2 - dx 2 - dy 2 . (39) 

(4M) 2 r v y I 

Parameter M plays for metric of the charged black hole the same role as the mass M for the 
neutral one. 

We now follow the arguments of ref. || for the Schwarzschild black hole. Define Rindler 
coordinates 

U = T + Z = pe tm y v = T _ z = _ pe -t/4M, (40 ) 



Then metric ( |39[ ) is simply Minkowski space in these coordinates 

ds 2 = dudv - dx 2 - dy 2 . (41) 

Let kn,, be the momenta of in-going and atmosphere (Hawking) particles respectively. The 
gravitational field of massless point-like particle in Minkowski space is described by the line 
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ds 2 = du 



x 2 



dv + 2k v In I — ^ I 5{u — uo)du 



dx 2 - dy 2 , (42) 



where x 2 = x 2 + y 2 , u = T + z and v = T — z. The massless particle moves in the v direction 
with constant uq and momentum k v . 
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The effect of the shock wave on the massless particle propagating in the metric (|42|) with 
initial momentum p^ is a discontinuity in the v direction at u = uq: 



x 2 



Av = - 2kVln [w)> (43) 

and a refraction in the transverse direction: 

Px(u) ~ Px = -^-Xp v 8(u - U ), (44) 

and similarly for p y (u) . 

By forming wave packets describing a high angular momentum Hawking particle before and 
after crossing the shock wave, and then calculating their scalar product one finds that the 
probability to be in the same state after crossing the shock wave is|| 

M 2 e 2 

P, 1 

1 p 4 ' 

where e is the energy of the atmosphere particle. Also, the probability for one particle in the 
atmosphere to have changed angular momentum is 

M 2 e 2 

Pai^o = 

This is the result of interaction of in-falling particle with one atmosphere particle. 

The number of particles which are affected by the shock wave of the in-going particle when 



it reaches p can be deduced from eqs. (19,21): 



N(p) ~ S(p) ~ -t 



Thus the probability for the atmosphere above p to remain in the initial state is 



/ tC/t2 2 \ N(p) 

Pf n f = R N{P) = [ 1 - — \ ~ r-M^rl/f 



tot ~ - ^ - —j£ 

We see that the state of the atmosphere is changed when the particle reaches p = (Mr + e) 1//3 . 

The minimal e one can consider is — since otherwise the wavelength of the in-going particle 
is larger than the radius of the black hole. Thus the minimal p is given by the following eq. 

/ Jl \ 1/3 



pr, 



which coincides with eq. (26). 

It can be easily shown using the same arguments as in ref. Q that the information carried by 
an in-going massless spin-less charged particle is encoded in the state of the atmosphere when 
the particle reaches p m i n . 
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4 Discussion 



We saw in the previous sections that the high angular momentum Hawking particles get reflected 
by the centrifugal barrier near the black hole horizon. It turns out that the number, total energy 
and other statistical quantities of these particles, which are found in the region over the horizon, 
depend only on the black hole parameters. This allows to speak about the atmosphere of 
Hawking particles near the horizon, as about pure quantum geometrical phenomenon. 

Since the Hawking radiation has a Planckian structure it is natural to say that a system 
consisting of a black hole and its radiation is in the thermal equilibrium state. One is interested 
to know thermodynamical parameters of the atmosphere in this state. In order to calculate 
entropy and the free energy of the atmosphere we counted the number of modes of the scalar 
field near the horizon, and then used well-known statistical thermodynamic formulae. Both 
entropy and the free energy are proved to be divergent at the horizon. 

In order to deal with these divergences we introduced a cutoff at the surface r = r + + h. The 
leading order term 0(/i -1 ) in the expression for the entropy is the contribution of the horizon to 
the total atmosphere entropy. The value of h was fixed by the requirement that the black hole 
entropy be a quarter of the horizon area, according to the Hawking-Bekenstein formula. This 
program was carried out for the Schwarzschild black hole in ref . , and for charged and rotating 



black holes in sections [2.1| and 2.3 respectively. We found that the value of the cutoff expressed 
in terms of the invariant distance is the same for all black holes independently of their mass, 
charge and angular momentum and is of order unity. At this (Planck) scale the semi-classical 
approach breaks down. 

Up to now the back reaction of the emerging radiation was neglected. Accounting for the 
back reaction gives rise to the black hole's mass decreasing and to the black hole's parameters 
fluctuation. The first phenomenon can be neglected at short time scales, provided that the black 
hole's energy is much larger than the Planck mass. The second one is responsible for the strong 
gravitational interactions occurring near the horizon. 

Indeed, by studying the atmosphere's mass fluctuations and interactions between Hawking 
and in-coming particles, we showed that the classical trajectories near the horizon cease to 
exist at invariant distances smaller than p m in = [ r +/( r + — r_)] 1//3 . This means that one cannot 
continue to use a non-quantized background metric at this scale, the semi-classical approach 
becomes invalid. 

If the black hole has small charge and angular momentum then p m in <C M and the question is: 
how does the high frequency part of the Hawking radiation spectrum affected by the gravitational 
interactions. Our analysis does not provide the analytical answer to this question. However, 
the fact that information is encoded out of the horizon, may be the reason pro the possibility of 
S"-matrix construction advocated by 't Hooft. 

An essentially new effect arises, when we turn our attention to the extreme black holes. As 
the distance between horizon r + and the surface r_ becomes smaller, the strong gravitational 
interactions occur at greater invariant distance, the smaller frequency Hawking quanta are af- 
fected, and so the whole spectrum becomes non-thermal. This in turn implies that in the case 
of the extreme black hole the information of an in-going particle is encoded at distances much 
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larger than black hole's radius. In this case, the semi-classical approach breaks down at all 
scales. 

Despite the fact that the minimum invariant distance from the horizon p m in at which the 
semi-classical theory still holds becomes infinite for the extreme black holes, its counterpart in 
the Boyer-Lindquist coordinates li mm = h(p m i n ) vanishes. However it still remains greater than 
value of the cutoff. Indeed, definition of the invariant distance implies that in general 

( r+ -r_)i/3 

ll>min ~ Pmin ~ M 2 / 3 

The cutoff in these coordinates is 

(r+ - r_) 
"-cutoff jr^ 

So, /icutoff < h m in for all black holes. Thus, one cannot reach the extreme black hole horizon 
(naked singularity) without entering a region below the h m i n which cannot be investigated with 
the semi-classical theory. We believe that the gravitational interactions which occur in this 
region somehow prohibit the particle to reach the horizon (and to violate causality). 

There is another simple argument which explains why one cannot reach the extreme black 
hole horizon in spite the fact that the non- invariant cutoff vanishes. The proper time it takes 
a freely falling observer to reach the horizon from some point R out of the horizon is infinite 
because this time is just the invariant distance 

M , 

At = p V2M - R, 

r + — r_ 

which is infinite for all R in our case. 

As we explained, the semi-classical approach is wrong if the radial coordinate in the Boyer- 
Lindquist coordinates is smaller than r + + h m i n . This leads to some inconsistency when we use 
the brick wall model, since we count semi-classical modes between r + + h cuto g and r + + h m i n 
also. However, we can neglect the contribution of this region, since the number of modes in the 
atmosphere is determined by the phase space volume in the nearest vicinity of the cutoff. T(p m i n ) 
is of order 0((M(r+ — r_)) 2 / 3 ) which is negligible compared to the leading order term 0(M 2 ). 
Moreover, the explicit (quantum gravitational) accounting for the processes occurring in this 
region should not change the estimate of p m in, since a completely different approach, — shock 
wave model, which does not use the counting of states, gives the same result. It is interesting to 
note that in the extreme black hole case the problematic region shrinks and has zero measure. 

We did not discuss so far how an external observer learns about the black hole's horizon 
fluctuations. We noted, that the black hole mass is fixed in his frame (apart of small decrease 
caused by the back reaction). So, in order to obtain information about the horizon fluctuations, 
the external observer should study trajectories of test particles which pass near the horizon. The 
fluctuating horizon may trap, with some probability, a particle with definite quantum numbers. 
And it may release another particle, which, in general, will not bear the same quantum numbers, 
since we have no information about the state of order or disorder of matter inside the black hole. 
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Therefore, the out-going geodesies which pass near the horizon will be thermally averaged and 
will differ from the classical geodesies. 

The external observer may also ask another observer who is found near the horizon, to send 
him information about the atmosphere's state. Then he must process this information bearing 
in mind that the Hawking temperature of the observer in the atmosphere is blue-shifted, and the 
light signals are red-shifted compared to those in his frame. It may be interesting to investigate 
this issue in more detail. 

The fact that the Hawking radiation is merely fluctuations of vacuum, may put in doubt our 
approach to the atmosphere as a system of real on-shell particles. As was pointed out in ref. @, 
it is expected to be a valid approximation if the S'-matrix ansatz of 't Hooft is correct. 

In summary, our analysis suggests that strong gravitational interactions occur near the black 
hole horizon at an invariant distance of order p m in = — f-)] 1 / 3 . At smaller distances 

the semi-classical approach breaks down. 



A Appendix. Shock wave on charge black hole background 



In the Appendix we shall show that the shock wave of a special form (48), generated by a 
point-like massless particle can propagate on the charged non-extreme black hole background. 

Recall the general result obtained by Dray and 't Hooft|12]: Given a solution of the vacuum 
Einstein equations of the form 

ds 2 = 2A(u,v)dudv + g(u,v)h ij (x' l )dx l dx j ; (45) 

if the following conditions hold 

A jV \ u=0 = = g >v \u=o, (46) 

^ Af _i^L f = 32 npA 2 S(x), (47) 
9 9 

where / = f(x l ) represents the shift in v, Af is the Laplacian of / with respect to the 2-metric 
hij and x is the transverse distance, then the shift in v at u = can be introduced so that 
the resulting space-time solves the field equation with a photon at the origin x = of the (x l ) 
2-surface and u = 0. The resulting metric is then described as follows: 

ds 2 = 2A(u, v + 6f)du(dv + efrdx*) + g(u, v + 6f)h ij dx l dx j . (48) 

Let us verify whether conditions of this statement are satisfied in the case we are interested 
in. At first, we introduce new coordinates u, v which are labels for outgoing and in-going, radial, 
null geodesies [@]. The geodesies are solutions of the following equation: 

ds 2 = = ^dt 2 - ^-dr 2 . (49) 
r z D 

Hence they are given by: 

u = t — r* , v = t + r*, (50) 
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where r* is defined by (|5|) and reads as follows 



r* = r + Mln^ + 2M2 - Q \n r -^±, (51) 
Dq (r + — r_) r — r_ 

and I?o is integration constant. 

The line element (|l]) in terms of new coordinates reads: 

ds 2 = ^dvdu - r 2 (d6 2 + sin 2 9dip 2 ). (52) 

It has obvious pathology at the horizon: D(r)\ r=r+ = 0. In order to remove it we change 
coordinates again Q. Using eqs (|5"0|),(|5"l|) one gets 

e («-fi)/(4M) = e r*/2M = e r/2M ^ (r - r_)2~ a , (53) 

where 

2M 2 - Q 2 



a 



2M(r+ -r_) 

and Dq is fixed using requirement that equation (^) must yield the correct result at Q = 0. 
Thus, new coordinates are defined as follows 

u= _ e -umM. v = e W*M 7 (54) 

where (3 is a certain constant to be fixed in such a way as to remove the pathology in the metric 
(H). We have 

n = M^ e . fe _ )M | +a)/3(r _ ^Mf-*)^ 
where eqs. (|^) and (|5l|) were used. Clearly, the value of /? we need is 

P = A) = Y~ — ' ( 55 ) 

2 +° 

Thus, by substitution this into ( [52] ) the line element reads 

ds 2 = ^Ml e -§KT( r - r-) x -^- a ^dudv - r 2 (d9 2 + sin 2 Odip 2 ). (56) 



In coordinates u,v,0, and ip metric of the charged black hole (56) has no pathologies. The only 
singularity at r = is a physical one. 

Functions A and g defined in eq. (45) are 



A= 1 -^e-^(r-r-)^-^°, (57) 
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g = -r . 

Note the useful relation following from eqs. (|53|) and fl54|): 



(58) 



2M 



r — r_ 



It shows that ^4 and 5 are functions of the product u ■ v, so 

.A(ti • v) >v oc u; g(u • v) >v oc it, 

and therefore conditions (f5j) are satisfied. 
Also, 

g u = e 2M (r — r_) v 2 ^°? 



r — r + 



2(r — M) — D(- + -^_) - (r - r_) + (r - r+)(a - -) + — 
V 7 r 2M ; V ' v +A 2 ; 2M 



Following ref. |12|] we arrange the coordinates so that the massless particle is at 9 = = u. 
Then, upon substitution of these relations into equation (f47|) one obtains 

Af - X{M,Q)f = -27tk(M,Q)5{9), (59) 

where 

for 2 . 



A(M,Q) = — 



2(r+ — M) — (r+ - r_) + 



2M 



k(M, Q) = 2 8 p^-(r + - r-) 1+ ( a -^°e-™ 
Po 

Note that A(M,0) = 1, and k(M,0) = 2 9 pe~ 1 M i . This is the result obtained by Dray and 't 
Hooft. 

We can find solution to eq. ( |59"| ) by expanding / in terms of spherical harmonics. However, 
since only spherical harmonics with m = contribute, it is suffice to expand / in terms of 
Legendre polynomials Pi(cos9). Then 

OO 7,1 

Since the asymptotic behavior of the Legendre polynomials at large I is 

P, ~ — 

this sum converges for all A such that < A < 1, i.e. for non-extreme black holes. A = 
corresponds to the extreme black holes \Q\ = M; in this case k(M,M) = 2 8 pM 3 and / diverges 
according to the following law: 

ext k(M,M) 1 



/ 



2X(M,M) y/M-\Q\' 
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